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We propose a way to accelerate adiabatic dynamics of wave functions in quantum 
mechanics to obtain a final adiabatic state except for the spatially uniform phase in 
any desired short time. We develop the previous theory of fast-forward (Masuda & 
Nakamura 2008) so as to derive a driving potential for the fast-forward of the adi- 
abatic dynamics. A typical example is the fast-forward of adiabatic transport of a 
wave function which is the ideal transport in the sense that a stationary wave func- 
tion is transported to an aimed position in any desired short time without leaving 
any disturbance at the final time of the fast-forward. As other important exam- 
ples we show accelerated manipulations of wave functions such as their splitting 
and squeezing. The theory is also applicable to macroscopic quantum mechanics 
described by the nonlinear Schrodinger equation. 

Keywords: atom manipulation, mechanical control of atoms, quantum 

transport 

1. Introduction 

The adiabatic process occurs when the external parameter of Hamiltonian of the 
system is adiabatically changed. Quantum adiabatic theorem (Born & Fock 1928; 
Kato 1950; Messiah 1962), which states that, if the system is initially in an eigen- 
state of the instantaneous Hamiltonian, it remains so during the process, has been 
studied in various contexts (Berry 1984; Aharonov & Anandan 1987; Samuel & 
Bhandari 1988; Shapere & Wilczek 1989; Nakamura & Rice 1994; Bouwmeester et 
al. 1996; Farhi et al. 2001; Roland & Cerf 2002; Sarandy & Lidar 2005; Du et al. 
2008). The rate of a change in the parameter of Hamiltonian is infinitesimal, so 
that it takes infinite time to achieve the final result in the process. 

In our previous paper (Masuda & Nakamura 2008), we investigated a way to 
accelerate quantum dynamics with use of an additional phase of a wave function 
(WF). We can accelerate a given quantum dynamics to obtain a target state in any 
desired short time. This kind of acceleration is called the fast-forward of quantum 
dynamics. 

The idea of the adiabatic process seems to be incompatible with that of fast- 
forward. But here we combine these two ideas, that is, we propose a theory to 
accelerate the adiabatic dynamics in quantum mechanics and obtain, in any desired 
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short time, the target state originally accessible after infinite time through the 
adiabatic dynamics. By using this theory we can find a driving potential to generate 
the target state exactly except for a spatially-uniform time-dependent phase like 
dynamical and adiabatic phases (Berry 1984). The fast-forward of the adiabatic 
dynamics makes an ideal transport of quantum states possible: a stationary wave 
packet is moved to an aimed position without leaving any disturbance at the end 
of transport. After the transport, the wave packet becomes stationary again, and 
is in the same energy level as the initial one. 

Before embarking upon the main part of the text, we briefly summarize the 
previous theory of the fast-forward of quantum dynamics (Masuda & Nakamura 
2008) . In Schrodinger equation with a given potential Vq = Vo (x, t) and nonlinearity 
constant c (appearing in macroscopic quantum dynamics) 

lh ^T = -7^V 2 * + ^o(x,i)*o-co|*o| 2 *o, (1.1) 
at 2mo 

^>o (x, t) is supposed to be a known function of space (x) and time (t) and is called 
a standard state. For any long time T called a standard final time, we choose 
&o(t = T) as a target state that we are going to generate. Let ^ a (x,t) be a 
virtually fast-forwarded state of ty (x,t) defined by 

|*a(*) >= |*o(at) >, (1-2) 

where a is a time-independent magnification factor of the fast-forward. The time- 
evolution of the WF is speeded up for a > 1 and slowed down for < a < 1 like a 
slow-motion. A rewind can occur for a < 0, and the WF pauses when a = 0. 

In general, the magnification factor can be time-dependent, a — a(t). The 
time-evolution of a WF is accelerated and decelerated when a(t) is increasing and 
decreasing, respectively. In this case, the virtually fast-forwarded state is defined 
as, 

|*a(*) >= |*o(A(t)) >, (1.3) 

where 

A(i) = / a(t')dt'. (1.4) 
Jo 

Since the generation of ^ a requires an anomalous mass reduction, we can not 
generate ^ a (Masuda & Nakamura 2008). But we can obtain the target state by 
considering a fast-forwarded state V&ff = *ff(x, t) which differs from ty a by an 
additional space-dependent phase, / = f(x,t), as 

* FF (i) = e l ^ a (t) = etf*o(A(t)). (1.5) 

The Schrodinger equation for ^ff is given by 

ih—^ = -o — V 2 *ff + Vff(x, t)^ FF - c |*ff| 2 *ff, (1.6) 
dt 2m 

where Vf f (x, t) is called a driving potential. If we appropriately tune the initial 
and final behaviours of the time dependence of a (the detail is shown later) , the 
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additional phase can vanish at the final time of the fast- forward Tp, and we can 
obtain the exact target state, that is, 

*ff(7f) = * (T), (1.7) 

where Tp is the final time of the fast-forward defined by 



= / « 
JO 



(t)dt. (1.8) 



(In the case of a constant a, Tf — T/a.) 

From equations (1.1), (1.4), (1.5) and (1.6) we obtain an equation for the addi- 
tional phase / (Masuda & Nakamura 2008) 

|* Q | 2 V 2 / + 2Rc[* Q V*;] • V/ 

+ (a-l)Im[* Q V 2 *;] =0, (1.9) 

and the driving potential of the fast-forward Vff 
V FF = „ % -»*-^(v/)« 

+Re[-(o - l)^-V 2 *„/*<, + 1— V/ ■ V*„/*J 
ztoo mo 

-(a-l)co|* Q | 2 . (1.10) 

With use of the phase 77 = r/(x, f) of the standard state \&o, / is given by (Masuda 
& Nakamura 2008) 

/(x,t) - (a(t)-l)T?(x,A(t)), (1.11) 

which satisfies equation (1.9) and determines Vff in equation (1.10). We impose 
the initial and final conditions for a(t): a must start from 1, increase for a while 
and decrease back to 1 with da/dt = at the final time of the fast-forward. Then 
the additional phase / vanishes in the initial and final time of the fast-forward. 
Once we have a standard state, we can obtain a target state in any desired short 
time by applying the driving potential with suitably tuned a(t). 

However, in the fast-forward of the adiabatic dynamics we shall use infinitely- 
large a. Then the expression of Vff in equation (1.10) and / in equation (1.11) 
should diverge. This difficulty will be overcome by regularization of the standard 
potential and states which will be described in Section 2. In Section 3 we apply 
the theory to the ideal transport by accelerating the adiabatic transport of a wave 
function in a moving confining potential. Ideal manipulations such as wave packet 
splitting and squeezing are shown as other examples of the fast-forward of adiabatic 
dynamics. Section 4 is devoted to summary and discussions. 



2. Regularization of standard state and driving potential for 

fast-forward 

(a) Difficulty involved in standard adiabatic dynamics 

We consider a dynamics of a WF, ty , under the potential: Vq = Vo(x,R(t)) 
which is adiabatically varies, where R = R(t) is a parameter in the potential which 
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is adiabatically changed from a constant Rq as 

R(t)=Ro+st. (2.1) 

The constant value e is the rate of adiabatic change of R(t) with respect to time 
and is infinitesimal, that is, 



dR(t) 



e, (2.2) 



dt 

e < 1. (2.3) 



The Hamiltonian is represented as 

H = ^- + V (x,R(t)), (2.4) 
and Schrodingcr equation for *S? is given as 

ih ^T = -^—^0 + ^o(x, i2(t))*o - c |*o| 2 *o, (2.5) 
dt 2mo 

where Co is a nonlinearity constant. If a system is in the n-th energy eigenstate 
at the initial time, the adiabatic theorem guarantees that, in the limit e — > 0, ^>o 
remains in the n-th energy eigenstate of the instantaneous Hamiltonian. Then *o 
is represented as 

* (x,t,i2(t)) = n (x,i?(t))e-*io tB nW*'))dt' e ir(t) ) (26) 

where £J„ = E n (R) and </>„ = <j> n (x., R) are the n-th energy eigenvalue and eigenstate 
corresponding to the parameter R, respectively, and T = F(t) is the adiabatic phase 
defined by 

r(t) = Vo y dxc ^«^"' (27) 

which is independent of space coordinates. 4> n fulfills 

■V 2 0„ + 1/ O (x,i?)0„ - co|^„| 2 </>„ =E n (R) ( f )n . (2.9) 



2m ( 



The second factor in the right hand side of equation (2.6) is called a dynamical 
phase factor, which is also space- independent. Such ideal adiabatic dynamics with 
e — > takes infinite time until we have an aimed adiabatic state (target state). 

Our aim is to accelerate an adiabatic dynamics <fi n (R(0)) — > <p n {R{T)) aside 
from the dynamical and adiabatic phases by applying the theory of fast-forward 
with infinitely-large magnification factor a, and obtain a target state 4> n (R(T)) in 
any desired short time, where T is a standard final time which is taken to be 0(1/ e). 
During the fast-forward, we will take the limit e^O, a^oo,T^oo and ae <~ 1. 
In this acceleration, we do not care about spatially uniform phase. (The spatially 
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uniform phase can be controlled by spatially uniform potential, if it is necessary.) In 
applying the fast-forward to the adiabatic dynamics, we have to choose a standard 
state and Hamiltonian. But there is some ambiguity in the choice of the standard 
state and Hamiltonian, because the state which we want to accelerate is represented 
in the limit e — > and we will fast-forward it with infinitely-large magnification 
factor. One might think that we can take \&o and Ho in equations (2.6) and (2.4) 
as a standard state and Hamiltonian. However, such idea is not adequate because 
the state in equation (2.6) is in an expression of WF in the limit e — > and does 
not satisfy Schrodinger equation up to 0(e) with small but finite e (Kato 1950; Wu 
& Yang 2005). In other words quantum dynamics in equation (1.1) with finite e 
inevitably induces nonadiabatic transition, but ^ m equation (2.6) ignores such 
transitions. Then the original theory of fast-forward (Masuda & Nakamura 2008) 
is not applicable as it stands. To overcome this difficulty, we shall regularize the 
standard state and Hamiltonian corresponding to the adiabatic dynamics. 

In the fast-forward of the adiabatic dynamics in the limit e — > 0, a — > oo and 
ae <~ 1, the standard final time T is chosen as T — O(^). In this case, a regularized 
standard state and Hamiltonian should fulfill the following conditions: 

1. A regularized standard Hamiltonian and state of the fast-forward should agree 
with Hq and ^>o except for space- independent phase, respectively, in the limit e — > 0; 

2. The regularized standard state should satisfy the Schrodinger equation corre- 
sponding to the regularized standard Hamiltonian up to 0(e) with finite e. 

Hereafter ^ e9 ^ and H^ 9 ^ are prescribed to a regularized standard state and 
Hamiltonian, respectively, which fulfill the conditions 1 and 2. 



(b) Regularized standard state 

Here we derive expressions for *Q res ' and H^ 9 ^ . Let us consider a regularized 
Hamiltonian H^ e9 ^: 

H (re g) = + v {re 3 ) ^ ^ 

Zmo 

where the potential V r < ' re9 ' 1 in the regularized Hamiltonian, is given as 

V {re9) (^t) = Vo(*,R(t))+eV(*,t). (2.11) 

V is a real function of x and t to be determined a posteriori, which is introduced to 

incorporate the effect of non-adiabatic transitions. It is obvious that H^ 9 ^ agrees 
with Ho in the limit e — > 0, that is, 

\imHlf ea \ X ,t) = H (x,R(t)). (2.12) 

Suppose the regularized standard state to be represented as 

y { re9) (x,t,R(t)) = cf> n ( X ,R(t))e-ifoEUR(t'))dt' ei ee(^t)^ (2 13) 



where 9 is a real function of x and t. <j) n and E n are the same as used in equation 

)■' 
0. 



(2.6). ~fy Q rc9 * > coincides with ^ except for the spatially uniform phase in the limit 
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Thus ^Q re9 ^ in equation (2.13) and Hjf e9 ^ in equation (2.10) fulfill the condition 
1 given in the previous subsection. 9 and V are chosen so that ^ res ' satisfies 
Schrodinger equation corresponding to H^ e9 ^ up to 0(e) to fulfill the condition 2. 

Here we shall obtain 9 and V. Schrodinger equation for ^f^ e9 ^ is represented as 

= - J^y 2 *^ + V^ re9 \ X} R(t))¥ re9) - c |^ re9) | 2 ^ re9) . (2.14) 
dt 2m,Q 

Substituting equations (2.11) and (2.13) into equation (2.14), we find 

dd> d9 h 2 

iKjge + Enfa-K—efa = - — [V 2 ^ + 2ieV 9 ■ V cfi n - e 2 {V 9) 2 (fi n 
OR dt 2m 

+ie(X7 2 9)<f> n ] + V4> n + eV4> n 

-cM 2 cf> n . (2.15) 

With use of equation (2.9) in equation (2.15), we have 

^TTET - = -^—[2iV9 • V0„ + i(V 2 9)(fi n ] + V<fi n , (2.16) 

OH at Zmo 

where we eliminated the term of the second order in e. Multiplying both sides of 
equation (2.16) by we have 

~ = i [2W - ' ve + Mv2e] + 4^ 2 - (2 - 17) 

Equation (2.17) is equivalently represented by decomposing it into real and imagi- 
nary parts as 

\<fi n \ 2 V 2 9 + 2Re[0„VC] • W + ^°Re[0 n ^] = 0, (2.18) 

—Im[(fi n V(fi n } ■ ^9 + + Im^^A,] + ^r\€\ = 0. (2.19) 

mo n OR dt 

9(~x.,t) should satisfy equation (2.18). Equation (2.18) can be derived also from the 
continuity equation for ty^ 9 ^ (See Appendix A). V is then given in terms of 9 and 
(fi n , that is, 

V = -H§ - Hm[^MJ - -Im[^] ■ W. (2.20) 

dt OR m (fi n 

Equation (2.18) indicates that 9 is a function which is not dependent on t explicitly, 
leading to 

d9__dRd9__ d9_ 

dt ~ dt OR ~ e dR~ 1 ■ ' 

Therefore, in our approximation to suppress the terms of 0(e 2 ), equation (2.20) is 
reduced to 

V = -hlm[^/cfi n } - — Im[^] • V6. (2.22) 

OR m (fin 

Thus we have obtained equations which 9 and V should satisfy so that the conditions 
of regularized standard state and Hamiltonian arc fulfilled. 
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(c) Driving potential for fast-forward 

Here we obtain the driving potential Vff to fast-forward the regularized stan- 
dard state 5'g re9 ' 1 . Such limit are taken as e — > 0, a — > oo and ae <~ 1 in the 
fast-forward. The final time of the standard dynamics T is taken to be 0(1/ e). The 
driving potential Vff = Vff(x, t) in equation (1.10) is explicitly represented as 

VffM = a y (^)( x ,i?(A(t)))-ft|(x,i)-^-(V.f) 2 

+z-^V/ ■ V*< refl) /*|> reff) ] - (a - l)coK re9) | 2 , (2.23) 
mo 

where / = /(x, i) is the additional phase of the fast-forwarded state f 

*ff(x,*) - e j/ ( x '*)< e9) (x,A(t)) (2.24) 

with A(i) defined by equation (1.4). 

As described in Section 1, the additional phase / given by equation (1.11) is 
not convenient for the fast-forward of adiabatic dynamics, because Vff in equation 
(2.23) would diverge due to the infinitely-large a. Instead we can take V/ as 

V/ = (a - l)eV0, (2.25) 

which also satisfies equation (1.9) (see Appendix B). So we can determine / such 
that 

/(x,t) = (a(t)-l)e0(x,A(t)). (2.26) 

With use of / in equation (2.26), we can avoid the divergence of the driving po- 
tential. By using equations (2.9), (2.11), (2.13) and (2.26) in equation (2.23), the 
driving potential Vff is written as 

h 2 

Vff = aeV + V - (a-1)- — e 2 (W) 2 

2mo 

t, da a kf i\ dd 
— n——e0 — n{a — lie— 

dt K ' dt 

/r 'a-l) 2 e 2 (V8) 2 + (a-l)E n (R(A)). (2.27) 



2m 



Noting — a (t"j £ d<L ; j n the limit e — > 0, a — > oo, and ea ~ 1 we can 

suppress the terms of 0(a p e q ) with q > p > in equation (2.27). While the last 
term in equation (2.27) would diverge, we omit it because it concerns only with 
space- independent phase of $j?f and has no effect on the dynamics governed by 
Vff- Then the driving potential reduces to 

V FF = aeV + V - H^e6 - ha 2 e 2 ^- - — a 2 e 2 (W6) 2 . (2.28) 
dt on 2m 
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From equation (2.28) the driving potential coincides with the standard potential 
Vo when as = and da/dt = 0. Thus, as should grow from and come back to at 
t = Tp, where Tp is the final time of the fast-forward defined by equation (1.8). On 
the other hand, as — 0(1) during the fast-forward. The change in the parameter 
R in equation (2.1), during t = and Tp, is given by 

AR := R(A(T F )) - R(0) = s [ a{t)dt. (2.29) 

Jo 

What we have to do for the fast-forward of an adiabatic dynamics is to obtain 9 
from a given adiabatic dynamics <fi n (R(t)) and then apply the driving potential Vpp 
given in terms of and <p n in equation (2.28). 



3. Examples 

We show some examples of the fast-forward of adiabatic processes in the linear 
regime (c = 0), by numerically iterating equation (1.6) with Vff(x,() in equation 
(2.28). (Examples in the nonlinear regime (cq ^ 0) will be reported elsewhere.) 
We can obtain the target state in the adiabatic process in any desired short time 
by applying the theory in Section 2. In the following examples, the magnification 
factor is commonly chosen for < t < Tp as 

2ir 

a(t)s = v cos(— * t + 7r) + v, (3.1) 

J-F 

where v is time average of a(t)s during the fast-forwarding, and the final time of the 
fast-forward Tp is related to the standard final time T as Tp = sT/v (see equation 
(1.8)). We take sT and Tp as finite, where s is infinitesimal and T is infinitely-large. 




o iy v 

t 7> 

Figure 1. Time dependence of a(t)e for < t < Tf- 

Namely we aim to obtain the target state in finite time, while the state is supposed 
to be obtained after infinitely long time T in the original adiabatic dynamics. The 
time dependence of a(t)s is shown in figure 1. 



Article submitted to Royal Society 



Fast-forward of adiabatic dynamics 



9 



(a) Fast-forward of adiabatic transport 

Suppose we have a stationary wave packet (WP) in a confining potential in 2 
dimensions and we want to transport the WP without leaving any disturbance on 
the wave function after the transport. If we move the confining potential rapidly 
to transport the WP fast, the WF is radically affected, and it will oscillate in the 
confining potential after the transport. However we can realize an ideal transport 
without leaving any disturbance at the end of the transport by applying the present 
theory of fast-forward of the adiabatic transport. Suppose there is a stationary 
WF, ip(x)e~^ Ent , in the confining potential U — U(x,y). Then let the potential 
adiabatically slide with infinitesimal velocity s in x-direction. The sliding potential 
V = V {x, y, t) is given by 



V = U(x-et,y). (3.2) 
The regularized WF in the potential is supposed to be given as 

< e9) - ^e-^-V^, (3.3) 



where 



K = i>(x-et,y)=ip(x-R(i),y), (3.4) 
R(t) = st. (3.5) 



From equation (3.4) we have 



d(j) n __ 0(j) n 
dR ~ dx ' 



(3.6) 



Substituting equation (3.6) into equation (2.18) or equation (A 5) we have 

h 



V-[|^| 2 V0] = 2Re[^„]. (3.7) 
mo ox 

It can be easily confirmed that 

= —x (3.8) 

fulfill equation (3.7) for any function <f) n . By substituting equations (3.6) and (3.8) 
into equation (2.22) wc find 

V = 0, (3.9) 

and the potential V$ re9 \x, y, t) in equation (2.11) in the regularized Hamiltonian 
H^ e9 ^ is given as 

Vt 9) =V = U(x-et,y). (3.10) 

It should be noted that we can also derive equations (3.8) and (3.9) from the 
Galilean transformation of the coordinates and omitting phase of WF in 0(s 2 ) 
which docs not affect the fast-forward. 
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With use of equations (3.8) and (3.10) in equation (2.28), we have 

V FF (x, y, t) = U(x - eA(t),y) ~ ^m Q ex - ^Ve 2 . (3.11) 

Since we have derived 9 without giving any specific form of 4> ni the formula of the 
driving potential in equation (3.11) is independent of the profile of the WF that we 
are going to transport. 

As an example of the fast-forward of the adiabatic transport, we choose a ground 
state in a harmonic potential: V g ' — {(x — xo(t)) 2 + (y — y ) 2 }. The centre 
of potential is adiabatically moved in x-direction as 

x (t) = x {0) + et. (3.12) 

xq corresponds to R in Section 2. The regularized standard state \E' res ' ) is 

*^(x,i)=^e^ ee (!^ ) lexp[-^{(x- a;o (t)) 2 + ( y - yo ) 2 } 

-iLut]e lee , (3.13) 

with 9 given by equation (3.8). We put sT = 20 and Tp — 6.25, where e is infinites- 
imal and T is infinitely-large. The centre of mass moves from (x,y) = (xo(0),yo) 
to (xq(0) + eT,yo) during the fast-forward, where yo is a constant. Figure 2 is the 




Figure 2. Spatio-temporal dependence of the driving potential (broken line) and the am- 
plitude of the fast-forwarded state |^ff| 2 (solid line) at y — yo. The parameters are taken 
as eT = 20, T F = 6.25, ^ = 1, u = 1, x (0) = 6 and y = 16. 

spatio-temporal dependence of the driving potential (broken line) and the ampli- 
tude of the fast-forwarded state (solid line) at y = yo- Figure 3 is spatio-temporal 
dependence of the real part of WF of the fast-forwarded state at y — yo- Due to the 
additional phase /, the real part of WF shows spatial oscillations during the fast- 
forward, but at the final time of the fast-forward such an oscillation disappears and 
$fjr agrees with the target state except for the spatially uniform phase. We can 
also fast-forward the adiabatic transport of WP in any excited state, and obtain the 
target state in finite time. This is the ideal transport of WPs because we can obtain 
the target state without any disturbances at the final time of the fast-forward. And 
WP becomes stationary again after the transport when ae = and % = 0. 
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Figure 3. Spatio-temporal dependence of the real part of wave function of the 
fast-forwarded state at y = yo- The parameters are the same as in figure 2. 

To numerically check the validity of the fast-forward, we calculate the fidelity 
defined as 

F = \<* FF (W£ e8 HMt))>\, (3-14) 

i.e. the overlap between the fast-forwarded state ^>FF{t) and the corresponding 
standard one ^ eg) (A(t)) , where ^ { Q re9) (A(7») = ^ { Q re9) (T). The time dependence 
of the fidelity defined by equation (3.14) is shown in figure 4. We find that the fidelity 
first decreases from unity due to the additional phase / of the fast-forwarded state, 
but at the final time it becomes unity with high precision (|1 — F\ < 1CP 5 ) again. 
We have thus obtained the adiabatically-accessible target state in a good accuracy 
in finite time Tf — 6.25. 




t T F =625 

Figure 4. The time dependence of the fidelity. 



(&) Fast-forward of wave packet splitting 

We then show the fast-forward of WP splitting in 1-dimension. Here we calculate 
9 by a numerical integration of equation (2.18) or equation (A 5). Let us consider 
the dynamics of a Gaussian wave packet (WP) to be split into a pair of separated 
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WPs. The regularized standard WF is chosen as 

<4 re9) = <i> n e le8 = h(R){(l- R)exp[-^x 2 } + i?x 2 exp[-^x 2 ]}e lee , (3.15) 

where R is defined by equation (2.1) with Rq = 0, and h(R) is a normalization 
constant represented as 

h(R) = [(1 -Rf^+il- R)R^- 2 + \r 2 ^]- 1/2 - (3.16) 

The potential in the regularized standard Hamiltonian is represented as 
v {reg) _ & Ra 2 x 4 + {(1 - R)a 2 - 5Ra}x 2 +2R-(1- R)a 

V ° ~ 2^~ Rx 2 + (i-R) ' ( } 

Note that V in equation (2.22) vanishes due to the absence of space-dependent 
phase of 0(1) (Im(^r/</>„) = Im(@±jr / 4> n ) = 0) in the regularized standard state 

in equation (3.15). It is easily confirmed that for constant R, "J^ 69 "* in equation 
(3.15) with 9 — is a zero-energy eigenstate under the potential V^ res ' ) in equation 
(3.17). While for R = the WF is a simple Gaussian, for R = 1 it becomes spatially- 
separated double Gaussians. The spatial distribution of |^ , [ ) re9 ' ) | 2 for various R from 
to 1 is shown in figure 5. The standard final time T is taken as T — 1/e and 




Figure 5. ^-dependence of ^o re9 ' | 2 f° r various R. 

R(0) = 0. As R is gradually changed from to 1, changes from a single-peaked 
WP at t = to a doubled peaked one at t = T, namely a splitting of a WP 
occurs. The final time of the fast-forward Tp is related to T as Tp = eT/v. 9 is 
obtained by numerical integration of equation (2.18) under the boundary conditions 
(a? = 0) = and 9{x = 0) = 0. The spatio-temporal dependence of 9 is shown in 
figure 6. 

The driving potential is calculated from equation (2.28) with use of equations 
(3.15) and (3.17) and 9. Vff (solid line) and I^ffI 2 (broken line) which is ac- 
celerated by Vff are shown in figure 7. The conversion from a hump to a hollow 
in the central region of the potential in figure 7 is caused by the deceleration of 
a(t)e which suppresses the splitting force. The fidelity defined by equation (3.14) 
is confirmed to be back to unity with a high numerical precision (0.999) at Tp (see 
figure 8). 
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Figure 6. Spatio-temporal dependence of 9. The parameters are chosen as eT — 1, 
T F = 0.025, = 1.0 and a = 1. 




Figure 7. Spatio-temporal dependence of |»T/ff| 2 x 10 5 (broken line) and Vff (solid 
line). The parameters are the same as in figure 6. 




t 7> = 0.025 

Figure 8. Time dependence of fidelity. 

(c) Fast-forward of wave packet squeezing: case ofV = 

The fast-forward of WP squeezing is also an important example of the accel- 
erated WF manipulation. Let us consider a WF with a Gaussian distribution in a 
harmonic potential in 1-dimension which is squeezed adiabatically. The regularized 
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standard state is represented as 



= fa™ ee (»V< exp[-^*»]e«. (3.18) 
Tin In 



The regularized standard potential is 

Vt 3) = (3.19) 
In the adiabatic dynamics, the potential curvature to is gradually increased as 

u(t) = lu + st, (3.20) 
where ujo is a constant. In this case, we can easily confirm that 

6 = -^ (3.21) 

satisfies equation (2.18), and that V in the regularized standard potential in equa- 
tion (2.22) vanishes, because lm(^j-/(f> n ) = Im(^p-/(/5 n ) — 0. With use of equations 
(3.18), (3.19) and (3.21) in equation (2.28) we obtain the driving potential as 

m uj 2 {A(t)) da mo Jmo 2 2 2 

- t 2 lTt £ Mm MAl ] ' ( } 

where A(t) is defined by equation (1.4). Spatio-temporal dependence of the driving 
potential Vff is shown in figure 9. The conversion of the potential curvature from 
positive to negative in figure 9 is caused by the deceleration of a(t)e which sup- 
presses the squeezing force. Time dependence of the fidelity is shown in figure 10. 
The fidelity becomes unity with a numerical precision of 0.999 at the final time of 
the fast-forward and the WP squeezing has been carried out in a short time. 

(d) Fast-forward of wave packet squeezing: case of V ^ 

In the previous examples, the right-hand side of equation (2.22) was always 
vanishing (for examples because of the absence of space-dependent phase of 0(1) 
in the standard WF). While in many practical cases we do not have to choose a 
stationary state with a space-dependent phase of 0(1) for a standard state, WF 
can in principle have such a phase. Let us finally show the fast-forward of WP- 
squeezing as an example of the fast-forward with non-zero V by having recourse 
to a stationary state with space-dependent phase 77 of O(l). Here the regularized 
standard state with E n = is chosen as 

= ^ ee (»V4 exp[-^x 2 ]e^e^, (3.23) 



with 77 given as 



V = B I cxp[^x' 2 }dx\ (3.24) 







where B is a real constant. The regularized standard potential that guarantees the 
eigenstate in equation (3.23) is written as 

Vt 3) = ~ + ^f* 2 ^ exp[^* 2 ] + eV. (3.25) 
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Figure 9. Spatio-temporal dependence of the driving potential (solid line) and the ampli- 
tude of the fast-forwarded state |^ff| 2 x 3.0 x 10 4 (broken line). The parameters are taken 
as w = 0.25, eT = 0.75, T F = 0.025 and = 1.0. The inset represents the potential 
(thick line) and wave function (thin line) distribution for —3 < x < 3 at t — 0. 




Figure 10. Time dependence of fidelity. 



lu is a parameter which is adiabatically changed as 

uj(t) =lu + et, (3.26) 

where ujq is the initial value of uj. ^^ e9 ^ with 6 = and constant ui in equation 
(3.23) stands for the zero-energy eigenstate trapped in the central hollow of the 
potential barrier V$ re9 ^ with V — (see the inset in figure 11). And the WP is 
squeezed adiabatically. 
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From equation (2.18), which is not affected by 77, 9 is obtained as 





m 2 
Afruj 



(3.27) 



Substituting equations (3.24) and (3.27) and <j) n in equation (3.23) into equation 
(2.22), V is obtained as 

V = -moB f x' 2 exp[^x' 2 ]dx' + A Bxe xp[^x 2 ]. (3.28) 
Jq ft 2lj h 

Substituting equations (3.27) and (3.28) into equation (2.28), we obtain the driving 
potential Vff, whose spatio-temporal dependence is shown in figure 11. The solid 



5x10 



-4x10 




Figure 11. Spatio-temporal dependence of Vff (solid lines). The broken lines represent the 
corresponding WF |*ff| 2 X 6.0 x 10 4 . The parameters are taken as ll> = 0.25, eT = 0.75, 
Tf = 0.025, = 1.0 and B — 0.1. The inset represents the distribution of the potential 
(thick line) and WF (thin line) at t = 0. 



and broken lines represent the driving potential and the corresponding amplitude 
of WF, respectively. The conversion of the potential curvature from positive to 
negative in figure 11 can be explained in the same way as in the case of the example 
(c). The amplitude of fast-forwarded WF is shown in figure 12. The inset in figure 
12 represents the time dependence of the fidelity for < t < Tp. We confirmed 
that the fidelity comes back to unity with a high precision at Tp- 



4. Conclusion 

We have shown the way to accelerate the adiabatic dynamics in microscopic and 
macroscopic quantum mechanics, with use of infinitely-large magnification factor 
of the fast-forward and regularization of standard states and Hamiltonian. One 
can obtain the target state exactly except for spatially-uniform (dynamical and 
adiabatic) phases in any desired short time, while in the adiabatic dynamics the 
target state is accessible after infinite time. A noble feature of the fast-forward of 
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adiabatic dynamics, which is distinct from our previous theory, is that we should 
regularize in advance the standard state so as to satisfy Schrodinger equation up to 
O(e) for a small but finite rate e of temporal change in Hamiltonian. To avoid the 
divergence of the driving potential, we used a new form of the additional phase / in 
equation (2.26) which differs from one in our previous work (Masuda & Nakamura 
2008). 

Typical examples of the fast-forward of adiabatic dynamics were shown, e.g., 
the wave function transport, and the wave packet splitting and squeezing. The fast- 
forward of adiabatic transport makes the ideal transport of wave functions possible. 
The wave packet is rapidly transported to the adiabatically accessible targeted 
position, with leaving neither disturbance nor oscillation after the transport: the 
wave function becomes stationary again in the confining potential at the end of the 
fast-forward. We confirmed that the fidelity first decreases from unity due to the 
phase / of the fast-forwarded state, but at the final time it becomes unity with high 
precision again. The examples of the wave packet squeezing and splitting show a 
way to fast-forward the adiabatic manipulation of wave packets without leaving any 
residual disturbance. The framework of this theory is applicable to the macroscopic 
quantum mechanics described by the nonlinear Schrodinger equation, and one can 
expect the high-speed manipulations of wave packets in Bose Einstein condensates, 
which will be reported elsewhere. 



S.M. thanks Japan Society for the Promotion of Science for its financial support. K.N. 
is grateful for the support received through a project of the Uzbek Academy of Sciences 
(FA-F2-084). We thanks D. Matraslov, B. Baizakov, B. Abdullaev, M. Musakhanov, S. 
Tanimura and S. Sawada for useful discussions. 



Article submitted to Royal Society 



18 



Shumpei Masuda & Katsuhiro Nakamura 



Appendix A. Derivation of equation (2.18) from continuity 

equation 

Equation (2.18) is also obtained from the continuity equation. Suppose that Schrodingcr 
equation for the regularized standard state \]/ re ^ is given by equation (2.14). And 
suppose ^ , p re9 ' ) to be represented by equation (2.13), where <p n satisfies equations 
(2.8) and (2.9). Then we have 

J t \< 9) \ 2 = d ^< 9) + ^ >d ^- = 2 £ Re[§§M (Al) 
where we used equation (2.8). On the other hand, we have the continuity equation: 

A|^)|2 = _i V . Im[^ re9) *V< es) ]. (A2) 
dt mo 

Substituting equation (2.13) into equation (A 2), we have 

||*<™»>| 2 = -A V .I m [^V0„ +J e|0„| 2 V0]. (A3) 
at too 

By the way we obtain 

V ■ Im^Wn] = ( A4 ) 

by multiplying 0* on the both sides of equation (2.9) and taking their imaginary 
part. Therefore by using equations (Al), (A3) and (A 4), we obtain 

— V-[|^ n | 2 Ve] = 2Re[^|0„]. (A5) 
too otx 

V# must satisfy this equation. From equation (A 5), we can reach equation (2.18). 

Appendix B. Additional phase of fast-forwarded state 

We show that the gradient of the additional phase / in equation (2.25) satisfies 
equation (1.9). , J , re9 ' in equation (2.13) includes the amplitude factor 

M^m) = \4>n\e iri ^ m) , (Bl) 

where r\ = ry(x, R{t)) is the space-dependent phase. With the use of equations (2.9) 
and (B 1) we have 

CV 2 0„ - 0„V 2 C = 2\<p n \[2Vr, ■ V|0„| + (V 2 r7)|^|] - 0. (B2) 
Thus, 7] satisfies 

2V|</> rl |Vr7+|0„|V 2 r/ = O. (B3) 

Noting in equation (1.9) is now given by \I/ a (t) = , J' re9 ' ) (A(i)) in the regularized 
case, we substitute equations (2.13) and (B 1) into equation (1.9), and obtain 

|0„| 2 V 2 / + 2|0„|V|0„|-V/ + (a-l)[-2|0„|V(»7 + e0)-V|0 n | 

-\^ n \ 2 V 2 (v + e9)} = 0. (B4) 

By using equation (B3) in equation (B4), we can easily check that V/ defined by 
equation (2.25) satisfies equation (B4), i.e., equation (1.9). 
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